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^ We consider the case of /{"-gravity and perform a detailed analysis of the dynamics 

, in Bianchi I cosmologies which exhibit local rotational symmetry (LRS). We find exact 

' solutions and study their behaviour and stability in terms of the values of the parameter 

n. In particular, we found a set of cosmic histories in which the universe is initially 
isotropic, then develops shear anisotropics which approaches a constant value. 
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1. Introduction 



The dynamical systems approach-*^ has been used with great success over the last 
. 30 years, to a gain (qualitative) description of the global dynamics of cosmologi- 

\ cal models. This method provides a useful tool for finding exact solutions which 

correspond to fixed points of the system. 

Carloni et aP have recently used this method to study the dynamics of i?"- 
theories in Friedmann-Lemaitre-Robertson- Walker (FLRW) universes. Clifton and 
Barrow^ used the dynamical systems approach to determine the extent to which 
exact solutions can be considered as attractors of spatially flat universes at late 
times. They compared the predictions of these results with a range of observations 
and argued that the parameter n in FLRW may only deviate from GR by a very 
rS ' small amount (n — 1 ^ 10~^^). 

' The main aim of this paper"* is to see how the shear behaves in LRS Bianchi I 

cosmologies in R"- gravity and whether these models isotropises at early and late 
times. To achieve this goal we use the theory of dynamical systems* to analyse the 
system of equations governing the evolution of this model with and without matter. 

2. Dynamics in LRS Bianchi I cosmologies 

The cosmological equations that we require for our analysis are: 

- 0, (1) 
= 0, (2) 

(3) 

A = -(1 + Hm0, (4) 





2a' 




-(n-l)|e- 








2n 


nR"-- 


1 o2 2 

36 -cr - 


\- (n - 












2n 






-f {n 








CT = - [ e 

















1 



February 6, 2008 17:8 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



2 



where 9 = ^ is the volume expansion and a is the shear (cr^ = \(Jah<^°'^\ 

In order to convert the above equations into a system of autonomous first order 
differential equations, we define the following set of expansion normalised variables; 

whose equations are 

n-l 

x' = 2/(2 + a;) ^—{2 + nx) -2x- 2x^ + xz + (1 - iw)z + 2xT,, (6) 

n ~ 1 

y' = — ^ [(3 - 2n)x ~ 2y + 2{n - \)z + 4(n - 1)E + 2(n - 1)1 , 
n — 1 



E' = 2 ( -2 + 2S ^ 2x + z ) S, 



2?7 

2z - (1 + 3w) - 3a; s!- + 4E 

n — 1 



z = z 
1 — S + a; — ?/ — z = 0, 

where primes denote derivatives with respect to a new time variable t = Ina. 
The solutions associated to the fixed points can be obtained from: 

e^^e^, «^(2 + E,-;^y,)"\ (7) 

- = -^e, /? = 2 + E, + y, + z,. (8) 

(7 3 

Under the condition that \ and the terms inside the brackets of ([71), do not add 
up to zero, these equations may be integrated to give the following solutions 

a = flQ — to)" , G = aoa^^ . (9) 

In terms of our expansion normalised variables (O, the energy density is given 

by 

fi (X zy''~^e^'\ (10) 

From this relation it can be seen that when z ~ and y ^ Q the energy density 
is zero. However when y — Q and z ^ the behaviour of does depend on the 
value of n. In this case the energy density is zero when n > 1 but is divergent when 
n < 1. When both y and z are equal to zero and n < 1, one can only determine the 
behaviour of /z by direct substitution into the cosmological equations. 



2.1. Vacuum analysis 

The vacuum case is characterised by z = 0. In this case we can obtain the fixed 
points of ([6]) by setting E' = and y' = 0. We find one isotropic fixed point 
(„4) : {q, l^l^f) ^"^^ ^ ^^^^ '^^ fixed points (£i) : (E*,0), with non-vanishing shear. 
The isotropic fixed point is an attractor (stable node) for values of the parameter 
n in the ranges n < 1/2, 1/2 < rt < 1 and n > 5/4. In the range 1 < n < 5/4 this 
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Fig. 1. Phase space of the vacuum model, (a) 1/2 < n < 1. (b) 1 < n < 5/4. The shaded regions 
represents the regions of initial conditions for which the shear will always evolve faster than in 
GR. 



point is a repeller (unstable node) and therefore may be seen as a past attractor. 
However, we also have attractors for (cr/i?)* << 1 on £i. Therefore inflation may 
not be needed since the shear anisotropy approaches a constant value which may 
be chosen as the expansion normahsed shear observed today^"^ ((cr/iJ)* < 10~^), 
provided that other observational constraints such as nucleosynthesis are satisfied. 

The phase space is divided into two regions by the line y = 1— S, which represents 
all points for which the shear dissipate at the same rate as in GR. The region 
y < 1 — T, represents a fast shear dissipation (FSD) regime where shear dissipates 
faster than in GR, and the region y > 1 — E, is a slow shear dissipation (SSD) 
regime where the shear dissipates slower than in GR (see Figure 1). 

2.2. Matter analysis 

Setting E' = 0, y' = and z' = wc obtain three isotropic fixed points and a line 
of fixed points with non- vanishing shear. 

We observe the same kind of behaviour as in the vacuum case; the phase space 
is however 3-dimensional, but is similarly divided into two regions, by the plane 
1 = T, + y + z. The space above the plane is a SSD region and below a FSD region. 

3. Conclusions 

In conclusion we have shown that R"- gravity modifies the dynamics of the shear in 
LRS Bianchi I cosmologies by altering the rate at which the shear dissipates. There 
are cases in which the shear always dissipates slower or faster than in GR, and there 
are ones which make the transition from first evolving faster and then slower (and 
vice versa) than in GR. 
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